UNIVERSALITY FOR ZEROS OF RANDOM ANALYTIC 

FUNCTIONS 



ZAKHAR KABLUCHKO AND DMITRY ZAPOROZHETS 

Abstract. Let £o>?ij--- be independent identically distributed (i.i.d.) ran- 
dom variables such that Elog(l + |£o|) < °°- We consider random analytic 
functions of the form 

oo 

Gn(z) = ^ £,kfk,nZ k , 
k = 

where fh,n are deterministic complex coefficients. Let v n be the random mea- 
sure assigning the same weight 1/n to each complex zero of G n . Assuming 
essentially that — — log/rt n i n ~^ U M as n — >■ oo, where u(t) is some function, 
we show that the measure v n converges weakly to some deterministic measure 
which is characterized in terms of the Legendre— Fenchel transform of u. The 
limiting measure is universal, that is it does not depend on the distribution 
of the £fc's. This result is applied to several ensembles of random analytic 
functions including the ensembles corresponding to the three two-dimensional 
geometries of constant curvature. As another application, we prove a random 
polynomial analogue of the circular law for random matrices. 



1. Introduction 

1.1. Statement of the problem. The simplest ensemble of random polynomials 
are the Kac polynomials defined as 

n 
k=Q 

where £oj£t . . . are non-degenerate independent identically distributed (i.i.d.) ran- 
dom variables. The distribution of zeros of Kac polynomials has been much stud- 
ied; see pEU [31 III [2S HI I2H1 EH H- :t is known that under a very mild mo- 
ment condition, the complex zeros of K„ cluster asymptotically near the unit circle 
T = {\z\ = 1} and that the distribution of zeros is asymptotically uniform with 
regard to the argument. Given an analytic function G which does not vanish iden- 
tically we consider a measure [Iq counting the complex zeros of G according to 
their multiplicities: 

[i G = ^2 n G (z)S(z). 

zeC:G(z)=0 
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Here, ng(z) is the multiplicity of the zero at z and S(z) is the unit point mass at 
z. Then, Ibragimov and Zaporozhets |14| proved that the following two conditions 
are equivalent: 

(1) With probability 1, the sequence of measures ^Mk„ converges as n — > oo 
weakly to the uniform probability distribution on T. 

(2) Elog(l + |£„|) < oo. 

Along with the Kac polynomials, many other remarkable ensembles of random 
polynomials (or, more generally, random analytic functions) appeared in the lit- 
erature. These ensembles are usually characterized by invariance properties with 
respect to certain groups of transformations and have the general form 

oo 
fe=0 

where £o,£i,.. . are i.i.d. random variables and fk, n are complex deterministic co- 
efficients. The aim of the present work is to study the distribution of zeros of G„ 
asymptotically as n — > oo. More precisely, we will show that under broad assump- 
tions on the coefficients fk, n , the random measure ^Mg„ converges, as n — > oo, to 
some limiting deterministic measure [i. The limiting measure [i is universal, that 
is it does not depend on the distribution of the random variables £kl see Figure [T] 
Universality has been much studied in the context of random matrices; see, e.g., 
|34j . The literature on random polynomials and random analytic functions usu- 
ally concentrates on the Gaussian case, since in this case explicit calculations are 
possible; see, e.g., [HI [3 HH EH ESI EH HH HI E01 [HI E]- The only ensemble of 
random polynomials for which universality is well-understood is the Kac ensemble; 
see [5H [U [TS1 [2] . In the context of random polynomials there have been a num- 
ber of results on the local universality in the distribution of zeros [T51 HH] . 
For example, universal character of local correlations between close zeros has been 
demonstrated for some models. In this work, our focus is different: we prove the 
universality of the distribution of zeros on the global scale. 

The paper is organized as follows. In Sections |2.1fj2.4| we state our results for 
a number of concrete ensembles of random analytic functions. These results are 



special cases of the general Theorem 2.9 whose statement, due to its technicality, 



is postponed to Section |2.5| Proofs are given in Sections [3j |4j [5j 

1.2. Notation. Let D r = {z E C : \z\ < r} be the open disk with radius r > 
centered at the origin. Let D = be the unit disk. Denote by A the Lebesgue 
measure on C. A Borel measure /i on a Polish space X is called locally finite (l.f.) 
if fJ,(A) < oo for every compact set A C X. A sequence \i n of l.f. measures on X 
converges vaguely to a l.f. measure \i if for every continuous, compactly supported 
function ip : X — >• K, 



(1) / ip(z)fj, n (dz) — > / cp(z)fx(dz). 

If fi n and n are probability measures, the vague convergence is equivalent to the 
more familiar weak convergence for which ([T]) is required to hold for all continuous, 
bounded <p. Let 9Jl(A) be the space of all l.f. measures on X endowed with the 
vague topology. A random measure on A is a random element with values in 
DJl(X). A sequence of random measures jj, n converges to a random measure /i in 
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Figure 1. Zeros of the Weyl random polynomial W„(z) = 
J2k=ot> k ^/k{ °f degree n = 2000. The zeros were divided by */n. 
Left: complex normal coefficients. Right: coefficients are real with 
P[log > t] = i -4 for t > 1. In both cases, the limiting distribu- 
tion of zeros is uniform on the unit disk. 



probability (respectively, a.s.), if ([!]) holds in probability (respectively, a.s.) for 
every continuous, compactly supported function ip. 

2. Statement of results 

2.1. The three invariant ensembles. Let ^0,^1, • ■ • be non-degenerate i.i.d. ran- 
dom variables. Fix a parameter a > 0. We start by considering the following three 
ensembles of random analytic functions (see, e.g., (3UD-2] anc * Figure [2]): 

[£Lo & ( "("- 1 )"-["- fc + 1 ) )%fc, (elliptic, n G N, z e C), 
F«W = j (^) a ^, (nat, n > 0,z € C), 

[ E » o& [ > n(n+i)...(n+*-l) ^° Jgfc) (hyperbolic, n > 0,z G D). 

Note that in the elliptic case F„ is a random polynomial of degree n, in the flat 
case it is a random entire function, whereas in the hyperbolic case it is a random 
analytic function defined on the unit disk D. 

In the particular case when a = 1/2 and are complex standard Gaussian 
with density z h-> 7r _1 exp{— |z| 2 } on C, the zero sets of these analytic functions 
possess remarkable invariance properties [31] [12]. Namely, in the flat case, the 
law of the zero set of F„ is invariant with respect to the rigid motions of C. In 
the elliptic (resp., hyperbolic) case, the law of the zero set of F n is invariant with 
respect to the isometries of the Riemann sphere C (resp., the unit disk D) preserving 
the spherical metric of constant positive curvature (resp., the hyperbolic metric of 
constant negative curvature). Thus, the three ensembles correspond to the three 
two-dimensional geometries of constant curvature. If a = 1/2 and are complex 
Gaussian, the expected number of zeros of F„ in a Borel set B can be computed 
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Figure 2. Zeros in the elliptic (left) and flat (right) models with 
a = 1/2 and n = 2000. The coefficients are complex normal. 

exactly: 

f = J B (1 + \z\ 2 )- 2 \{dz), (elliptic case, B C C), 
E[fi Fn (B)}= l^X(B), (flat case, B C C), 

[ J /_(1 - \z\ 2 )- 2 \(dz), (hyperbolic case, B C D). 

The next theorem states the universality for the distribution of zeros of F n . 

Theorem 2.1. Let £o>£lj--' ^ e non- degenerate i.i.d. random variables such that 
Elog(l + |£o|) < oo. As n — > oo, the sequence of random measures ^//f„ converges 
in probability to the deterministic measure having a density p a with respect to the 
Lebesgue measure, where 

{2^|z|i" 2 (l + |z|i)- 2 , [elliptic case, z € C), 
2^M"~ 2 , [flat case, z eC), 

2^M = ~ 2 (1 - \z\i )- 2 , [hyperbolic case, z € D). 

2.2. Littlewood Offbrd random polynomials. Next we consider an ensemble 
of random polynomials which was introduced by Littlewood and Offord 19, 2D]- It 
is related to the flat model. First we give some motivation. Let £o>£i>--> t> e non- 
degenerate i.i.d. random variables. Given a complex sequence wq,wi, . . . consider 
a random polynomial W„ defined by 

n 

(2) w n (z) = J2tkmz k - 

For Wk = 1 we recover the Kac polynomials, for which the zeros concentrate near 
the unit circle. The next result shows that the structure of the zeros does not differ 
essentially from the Kac case if the sequence Wk grows or decays not too fast. 

Theorem 2.2. Let £o>£l>--' be non- degenerate i.i.d. random variables such that 
Elog(l + |£o|) < oo. If lim.fc_j.oo rlog|u>fc| = w for some constant w £ R, then 
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Figure 3. Left: Zeros of the Szego polynomial s n (z) — ^2k=o T-i 
of degree n = 200. Right: Zeros of the Littlewood-Offord random 
polynomial W„(z) = ^fe=o £ fc tr °f degree n = 2000 with complex 
normal coefficients. The zeros were divided by n. 



the sequence of random measures ~/iw ra converges in probability to the uniform 
probability distribution on the circle of radius e~ w centered at the origin. 

We would like to construct examples where there is no concentration near a 
circle. Clearly, the sequence log|ii>fc| has to grow or decay superlinearly. We will 
consider decaying sequences, since for growing sequences the zeros concentrate near 
the origin. At a first attempt, it is natural to look at the case in which log \uik\ is 
a multiple of — fclogfc. More precisely, we make the following assumption on the 
sequence w^: 

(3) log \w k \ = — a(fclogfc — k) — (3k + o(k), k — > oo, 

where a > and ft G K are parameters. Particular cases are polynomials of the 
form 

w r W(z) = y -^-zK, w^(z) = y & **, wi 3 )(z) - y - ^ **. 

n v ; ^ (k\) a n w ^ k ak n w ^T(ak + l) 

k=o v ' fe=o fe=o v ; 

The family Wn has been studied by Littlewood and Offord [THl HO] m one of the 
earliest works on random polynomials. They were interested in the number of real 
zeros. In the next theorem we describe the limiting distribution of complex zeros 
ofW„. 

Theorem 2.3. Let £oj£ij-- - be non- degenerate i.i.d. random variables such that 
Elog(l + |£o|) < oo. Let u>o,wi,... be a complex sequence satisfying O). With 
probability 1, the sequence of random measures j/iw„ (e^n a -) converges to the de- 
terministic probability measure having the density 

(4) z ^ — 1 -|z|»- 2 l zeD 

lira 

with respect to the Lebesgue measure on C. 
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Remark 2.4. The measure nw n (e^n 01 -) counts points of the form / a , where z is 
a zero of W n . 

For the so-called Weyl random polynomials W„(z) = X)fe=o £ fc having a = 
1/2 and j3 = 0, the limiting distribution is uniform on D; see Figure [TJ This re- 
sult can be seen as an analogue of the famous circular law for the distribution of 
eigenvalues of the non-Hermitian random matrices with i.i.d. entries [33] [5]. For- 
rester and Honncr ffl[ stated the circular law for Weyl polynomials and discussed 
the differences and similarities between the matrix and the polynomial cases; see 
also [17]. The analogy between the random matrices with i.i.d. entries and the 
Weyl polynomials is not merely the coincidence of the limiting distributions. Both 
models are closely con necte d to the logarithmic potential theory on C with external 

for more details. 



2.6 



field ^\z\ 2 ; see Section 

Under a minor additional assumption on the coefficients Wk we can prove that 
the logarithmic moment condition is not only sufficient, but also necessary for the 
a.s. convergence of the empirical distribution of zeros. It is easy to check that the 
additional assumption is satisfied for W„ = w['' with £ = 1,2, 3. 

Theorem 2.5. Let£oi£i, ■ ■ ■ be non- degenerate i.i.d. random variables. LetWo,Wi, . . 
be a complex sequence satisfying ^ and such that for some C > 0, 

(5) \w n - k /w n \ < Ce pk n ak for allneN,k< n. 
Then, the following are equivalent: 

(1) With probability 1, the sequence of random measures -/Jw„ (e^n"-) con- 
verges to the probability measure on D with density Q. 

(2) Elog(l + |£„|) < oo. 

It should be stressed that in all our results we assume that the random variables 
are non-degenerate (that is, not a.s. constant). To see that this assumption is 
essential, consider the deterministic polynomials 

n 

(6) ^) = EF 

k=a 

A classical result of Szego [33] states that the zeros of s n (nz) cluster asymptotically 
(as n — > oo) along the curve {|ze 1_z | — 1} n D; see Figure [3] (left). This behavior 
is manifestly different from the distribution with density l/(27r|z|) on D we have 
obtained in Theorem [2^3] for the same polynomial with randomized coefficients; see 
Figure [3] (right). 

2.3. Littlewood Offord random entire function. Next we discuss a random 
entire function which also was introduced by Littlewood and Offord [21] [22] . Their 
aim was to describe the properties of a "typical" entire function of a given order 
1/a. Given a complex sequence Wo, Wi,... satisfying ^ consider a random entire 
function 

00 

(7) W(z) = Yj k w k z k . 



k=0 



Examples are given by 



w^hVttS-a w^( Z ) = y^z k , w< 3 >(--^ - ? 



^ (k\) a ' w ^ k ak ' w ^r(afc + l) 

fe=0 v ' fc=0 fc=0 v ; 
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The first function is essentially the flat model considered above, namely ~W(n a z) — 
F„(z). For a — 1 it is a randomized version of the Taylor series for the exponential. 
The last function is a randomized version of the Mittag-Leffler function. Our aim 
is to describe the density of zeros of W on the global scale. We have the following 
strengthening of the Hat case of Theorem |2.1| 



Theorem 2.6. Let £o>£ij--' be non- degenerate i.i.d. random variables such that 
Elog(l + |£o|) < oo. Let Wo, Wi, .. . be a complex sequence satisfying ([3]). With prob- 
ability 1, the random measure ^/Av converges to the deterministic measure 
having the density 

(8) Z ^277a N "~ 2 
with respect to the Lebesgue measure on C. 

As a corollary we obtain a law of large numbers for the number of zeros of W. 

Corollary 2.7. Let N(r) = /zw(D r ) be the number of zeros of W in the disk D r . 
Under the assumptions of Theorem \2. 



N(r) — e '<* r <* (1 + o(l)) a.s. asr-> oo. 



In the case a — 1/2 the limiting measure in Theorem 2.6 has constant density 
1/tt. The difference between the limiting densities in Theorem 2.3 and Theorem 2.6 
is that in the latter case there is no restriction to the unit circle. In the case of 
Bernoulli-distributed the function Q has been considered by Littlewood and 
Offord [211 [22]. Under the assumption log|ii)fc| ~ — aklogk as k — > oo they proved 
some estimates for the distribution of zeros of W. These estimates does not seem 
to imply Theorem |2.6[ since they are true up to multiplicative constants only. Also, 
the estimates of Littlewood and Offord are stated in terms of some more advanced 
functions counting the zeros and are not easily translated to our setting. 

Let us again stress the importance of the non-degeneracy assumption. The 
exponential function e z has no complex zeros, whereas the zeros of its randomized 
version X^o£ fc lT nave t ne global-scale density l/(27r|z|) on C. For the absolute 
values of the zeros, the limiting density is constant and equal to 1 on (0,oo). 

2.4. Randomized theta function. Given a parameter a € (0,1) U (l,oo) we 
consider a random analytic function 



££Lo&e nl " fe " ' z k , (case a < 1, z e D), 



Hn(2) 



Y,T=o^e n ak " zk -, (case a > l,z e 



Theorem 2.8. Let £o> £!.>••• ^ e non- degenerate i.i.d. random variables such that 
Elog(l + |£o|) < °°- As n — > oo, the sequence of random measures ^Mh„ converges 
in probability to the deterministic measure having the density 



1 1 

Z ^ 2na\l-a\ \zf 



logM 

a 

with respect to the Lebesgue measure on C. The density is restricted to D in the 
case a < 1 and to C\D in the case a > 1. 

As the parameter a crosses the value 1, the zeros of H„ jump from the unit disk 
D to its complement C\D; see Figure gj (left) for the case a = 2. Note that the 
case a — 1 corresponds formally to Kac polynomials for which the zeros are on 
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Figure 4. Left: Zeros of the randomized theta function H„(z) = 

k 2 

Sfc°=o£ fee ~~ zk with n = 2000 and normally distributed coeffi- 



cients. Right: Zeros of the random polynomial from Example 2.11 
with n — 300 and Cauchy-distributed coefficients. 



the boundary of D. The special case a — 2 corresponds to the randomized theta 
function 



DO 



(9) H n (z) = ]Tae-^z fe . 

fe=0 

The limiting distribution of zeros has the density g^fera on C\D. One can also take 
the sum in ^ over k € Z in which case the zeros fill the whole complex plane with 
the same density. 

A similar model, namely the polynomials Q„(z) = Y^k=o £fc e_fe z > where a > 1, 
has been considered by Schehr and Majumdar [35] • Assuming that ^ are real- 
valued they showed that almost all zeros of Q„ become real if a > 2. In our model, 
the distribution of the arguments of the zeros remains uniform for every a. 

2.5. The general result. We are going to state a theorem which contains all 
examples considered above as special cases. Let £o)£i> • • ■ be non-degenerate i.i.d. 
random variables such that Elog(l + |£o|) < oo. Consider a random Taylor series 

oo 

(10) G n (z)=J2^kfk,nZ k , 

fe=0 

where fk, n G C are deterministic coefficients. We assume that there is a function 
/ : [0, oo) — > [0, oo) and a number T £ (0, oo] such that 

(Al) f(t) > for t < T and f(t) = for t > T . 

(A2) / is continuous on [0, To), and, in the case To < +oo, left continuous at To. 

(A3) lim^oo sup ke[0An] \ \f k ,n\ 1/n - /(* )| = for every A > 0. 

(A4) R := liminf^^ fit)- 1 / 1 e (0, oo] and liming ^ \f k , n \~ 1/k > Ro- 



It will be shown later that Condition (A4) ensures that the series ( 10 ) defining 



G„ converges with probability 1 on the disk D flo . Let I : K — > K U {+00} be the 
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Legendre-Fenchel transform of the function u(t) = — log /(f), where logO = — oo. 
That is, 

(11) J(s) = sup(sf - u{t)) = sup(st + log f(t)). 

t>o t>o 

Note that / is a convex function, I(s) is finite for s < logi?o and J(s) = +00 for 
s > logi?o- Recall that ^g„ is the measure assigning to each zero of G„ a weight 
equal to its multiplicity. 

Theorem 2.9. Under the above assumptions, the sequence of random measures 
~HGn converges in probability to some deterministic locally finite measure fi on the 
disk V>r ■ The measure /i is rotationally invariant and is characterized by 

(12) n{D r ) =J'(logr), r<R . 

By convention, I' is the left derivative of I. Since I is convex, the left derivative 
exists everywhere on (—00, log Rq) and is a non-decreasing, left-continuous function. 



Since the supremum in ( 11 ) is taken over t > 0, we have lim, s _j._ 00 I'(s) = 0. Hence, 
fj, has no atom at zero. If I' is absolutely continuous on some interval (logri, logr2), 
then the density of fi on the annulus r% < \z\ < r 2 with respect to the Lebesgue 
measure on C is 

(13) M - «- 

2tt\z\ z 

It is possible to give a characterization of the measure \x without referring to 
the Lcgcndrc- Fcnchcl transform. The radial part of /1 is a measure /2 on (0, 00) 
defined by /i((0, r)) = /i(B r ). Suppose first that u is convex on (0, To) (which is the 
case in all our examples). Then, ft is the image of the Lebesgue measure on (0, 00) 
under the mapping t i-> e u where u' is the left derivative of u. This follows 
from the fact that (u') <_ = /' and (7')' <_ = u' by the Legendre-Fenchel duality, 
where <p*~(t) = mi{s E M : ip(s) > t} is the generalized left-continuous inverse of 
a non-decreasing function (p. In particular, the support of fi is contained in the 
annulus 

| e limt4.o«'(t) < | z | < e lim tt T «'(t)| 

and is equal to this annulus if v! has no jumps. In general, any jump of v! (or, 
by duality, any constancy interval of I') corresponds to a missing annulus in the 
support of fi. Also, any jump of I' (or, by duality, any constancy interval of u') 
corresponds to a circle with positive /z-measure. More precisely, if I' has a jump 
at s (or, by duality, u' takes the value s on an interval of positive length), then /i 
assigns a positive weight (equal to the size of the jump) to the circle of radius e s 
centered at the origin. In the case when u is non-convex we can apply the same 
considerations after replacing u by its convex hull. 

One may ask what measures \x may appear as limits in Theorem |2.9| Clearly, \x 
has to be rotationally invariant, with no atom at 0. The next theorem shows that 
there are no further essential restrictions. 

Theorem 2.10. Let fi be a rotationally invariant measure on C such that 

(1) fi{C\D Ro ) = 0, where R := sup{r > : /x(B r ) < 00} <E (0,oo]. 

(2) Jq ' fi(D r )r~ 1 dr < 00 for some (hence, every) R < R a . 



Then, there is a random Taylor series G n of the form { 10 ) with convergence radius 



a.s. Rq such that ^AtG„ converges in probability to [i on the disk Dr . 
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Example 2.11. Consider a random polynomial 

n 2ri , / _ s n 3n , 

(14) GbW =5:^ + 2» e *+(§) e Mf) • 

fc=0 fc=n+l ^ ' fc=2n+l 

We can apply Theorem |2.9| with 

'0, *e[0,l], (0, s<0, 

(log2)(t-l), fe[l,2], I()= \«, sG[0,log2], 

S (log3)i-log§, iG[2,3], U I 2s -log 2, s G [log 2, log 3] , 

+oo, f > 3, [3s — log 6, £>log3. 

The function u' has three constancy intervals of length 1 where it takes values 
0, log 2, log 3. Dually, the function I' has three jumps of size 1 at 0, log 2, log 3 and 
is locally constant outside these points. It follows that the limiting distribution of 
the zeros of G n is the sum of uniform probability distributions on three concentric 
circles with radii 1,2,3; see Figure [4] 

Remark 2.12. Suppose that G„ satisfies the assumptions of Theorem |2.9| Then, 
so does the derivative G^ (and, moreover, / is the same in both cases). Thus, the 
derivative of any fixed order of G„ has the same limiting distribution of zeros as G n 
itself. Similarly, for every complex sequence c„ such that lim^oo ^ log |c„| < /(0), 
the function G„(z) — c n satisfies the assumptions. Hence, the limiting distribution 
of the solutions of the equation G„(z) = c„ is the same as for the zeros of G„. 

2.6. Connection with logarithmic potentials and orthogonal polynomials. 

The identity ^Alog \z — Zo\ = 5{zq) implies the crucial formula 

(15) MGn = ^Alog|G„( 2 )|, 

where the Laplacian A should be understood in the distributional sense. The first 
step in the proof of Theorem |2.9| is to compute the limiting logarithmic potential 
in (pi 



Theorem 2.13. Under the assumptions of Section 2.5, for every z € Dfl o \{0}, 
(16) p n (z) := - log |G„(z)| A /(log \z\). 

71 n— foo 



Remark 2.14. Taking (151 and (16) together we obtain, at least formally, the fol- 
lowing equivalent version of (12): fx = ^A(/(log \z\)). 

The measure [i can be interpreted as an equilibrium measure for the logarithmic 
potential in the presence of an external field; see [24] for an account of logarithmic 
potential theory. Assume that the assumptions of Section [2~5| hold and additionally, 
the function u is convex on (0, To). Since we would like to deal with finite measures, 
let us take some k £ (0,To) and consider a truncated version of G n : 

[kti] 

then so does G^ with 



If G n satisfies assumptions (A1)-(A4) of Section 



2.5 



u(t), te[0,K] 

+00, t > K. 
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By Theorem 2.9 the measure ^M g m converges in probability to a limiting measure 
denoted by /J- K K 

Proposition 2.15. The measure /j,^ is the unique minimizer of the functional 

= 5 / / log , v{dz)v{dw) + I I(log\z\)v{dz) 
ZJcJc \z-w\ Jc 

in the set of all measures v on C which have total mass k. 

Denote by the support of the measure fx^ K \ From the interpretation of 
as the image of the Lebesgue measure under the mapping t i-> exp{(u( K ) )'(£)} the 
following monotonicity property follows: if Ki < K2, then S^ 1 ^ C S^K If the 
left derivative of has no constancy intervals, then the measures jjS Kl ^ and /i^ 2 ^ 
coincide on S^ Kl \ Thus, as k grows, the zeros fill larger and larger domains without 
changing their density on domains which are already occupied by the zeros. 

Let us finally mention a connection between the random polynomial models con- 
sidered above and orthogonal polynomials. Let to be a finite measure on C having 
an infinite support and finite moments of all orders. We can construct a sequence 
of polynomials p (z),pi(z), . . . which form an orthonormal system in the Hilbert 
space L 2 (m) and such that the degree of pk is k. At least for some special classes of 
measures to the empirical distribution of zeros \[ip h of the deterministic polynomial 
Pk converges, as k — > 00, to some limiting measure which can be interpreted as the 
equilibrium measure (for the logarithmic potential without external field) on the 
support of m. The simplest example are Chebyshev (or, more generally, Jacobi) 
polynomials whose zeros have limiting arcsine distribution. However, the equilib- 
rium property of the limiting empirical distribution may fail even for very simple 
measures to. For example, if to is rotationally invariant, then pk is a multiple of z k 
having a unique zero of multiplicity k at 0. 

We are going to show how the potential theoretic interpretation of the limit- 
ing distribution of zeros can be restored if we pass to random polynomials. Let 
wiii m 2, •■■ be finite, rotationally invariant measures on C having moments of all 
orders. The polynomials which are orthonormal with respect to to„ have the form 
Pk,n(z) = fk : ? l z k , k g N , where 

-1/2 

(17) h, n = ( / \z\ 2k m n (dz) 



Consider a random linear combination of these polynomials 

[kti] 



k,nZ 



k=0 

Some of the ensembles of random polynomials we considered above are particular 
cases of this setting. For example, taking m n (dz) — ^e~™' 2 ' \(dz) we have pk, n — 
(zy/n) k j\Jk\ thus recovering the Weyl polynomials. Kac polynomials are recovered 
if we take m n — At, the uniform probability distribution on T. See Table [l] for 
more examples. Observe that in all examples listed in Table [T] the sequence m n 
satisfies a large deviation principle (see [6]) with speed n and a certain rotationally 
invariant rate function 2Q(|z|). This means that 

limsup - logm n (K) < -2Q(K), liminf - logmjt/) > -2Q(U) 

n^oo n n->oo n 
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Table 1. Random analytic functions associated to orthogonal 
polynomials with rotationally invariant weights 



Model 


Pk,n 






Q(W) 


„("> 




Kac 


z fc 




At 


Jo, z 6 T, 
\+oo, z £ X 




T 


Flat (Weyl) 






5 ,-»w a 




1\ 


S ^ 
















Elliptic 


/t»(n-I)...(n-fc4 




1 n+l 


llog(l+ |z| 2 ) 






V fei 




- (l+| s |2)n+2 


- (l + |z| 2 ) 2 


















Hyperbolic 


/n(»+l)...(»+k- 




^(l-l-l 2 ) — 2 


-ilog(l - |Z| 2 ) 


1 1 




(z e d) 


V fei 




' (l-l-l 2 ) 2 


6 r ^=— 


Theta 

(z g C\D) 


e 2n z fc 




7r~3^i| z |— ™ lo 6 W 


(l°gM) 2 


1 





for every closed set K c C and every open set U C C. Here, Q(^4) = inf 2gj 4 Q(|z|) 
for any set A C C. For example, in the case of Weyl polynomials we have Q(|z|) = 
||z| 2 . On the other hand, in the case of Weyl polynomials, the external field in 
Proposition 2.15 is given by /(log |z|) = h\z\ 2 . 



This coincidence is a particular case 
of a more general statement, as we will argue now. We give the main idea skipping 



technical details. By Varadhan's lemma applied to (17 1 we have for all t > 0, 
1 



(18) lim - log/ tn ,„ = sup(£logr - Q(r)) = sup(ts - Q(e s )) =: u(t). 

n^oo n r >0 sER 

Imposing add itional minor a ssump tions it is possible to check that we can apply 



2.9 



and Proposition 



2.15 



to G ( n K) 



with I being the Legendre-Fenchel dual 

We have thus 
Some special 



Theorem 

of u. However, the Legendre-Fenchel dual of u is Q(e s ) by (18) 
identified the external field /(log|z|) with the rate function Q(\z\ 
cases are listed in Table [T] 

In the case of a not necessarily rotationally invariant weight m n = m which does 
not depend on n, the zeros of random combinations of orthogonal polynomials (with 
Gaussian coefficients) have been studied in [2H1 H] • Under appropriate conditions 
on m it has been shown that the zeros of such polynomials are asymptotically 
distributed according to the equilibrium measure on the support S of m. The 
external potential in this case is equal to +oo outside S and is on S. The external 
potential just restricts the equilibrium measure to S. This setting includes the Kac 
polynomials but it does not include other examples considered here. 

2.7. Open questions. We established the global universality for the distribution 
of complex zeros under a logarithmic moment condition. We strongly believe that 
local universality for complex zeros, as well as local and global universality for real 
zeros for the models considered here can be proved under stronger moment con- 
ditions. For example, there should be local universality for distributions from the 
domain of attraction of an a-stable law, a € (0,2], however, different values of a 
should correspond to different local universality classes. The distribution of zeros of 
the Kac polynomials K n in the case of logarithmically decaying tails has been stud- 
ied in [TB]. We believe that in the general setting considered here, similar results 
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should hold. In particular, the zeros should have a circle structure similar to that 
found in [15]; see Figure [l] (right). It should be possible to generalize our results to 



(systems of) random analytic functions in many variables. In Theorems 2.3 2.6 we 
proved the a.s. convergence for the empirical distribution of zeros in Littlewood- 
Offord models. We used the natural nested structure of these models. It is open 
whether convergence in probability can be replaced by the a.s. convergence in The- 
orem |2.1| We don't know whether the approach of Section |2.6| can be extended 
to polynomials which are orthogonal with respect to not necessarily rotationally 
invariant weights. 



3. Proofs: Special cases 
We are going to prove the results of Section [T] We will verify the assumptions 



of Section 2.5 and apply Theorem 2.9 Recall the notation u(t) — — log f(t) 



Proof of Theorem \2.2\ We can assume that w = since otherwise we can consider 
the polynomial W n (e~ w z). It follows from lirm-^oo rlog|ti>fc| = that assump- 



tions (A1)-(A4) of Section 2.5 are fulfilled with T = 1, Rq = +oo and 



1, *€[0,1], 0, ^[O,!], 

w [0, t>l, w \+oo, t>l. 

The Legendre-Fenchel transform of u is given by I[r) = max(0,r). It follows 



from (12 1 that [i is the uniform probability measure on T. □ 



Remark 3.1. Under a slightly more restrictive assumption Elog|£o| < oo, Theo- 
rem 2.2 can be deduced from the result of Hughes and Nikeghbali [T3] (which is 
partially based on the Erdos-Turan inequality). This method, however, requires a 
subexponential growth of the coefficients and therefore fails in all other examples 
we consider here. 



Proof of Theorem 2.1 By the Stirling formula, logn! = n log n—n+o(n) as n — > oo. 



It follows that assumption (A3) holds with 

{a(tlogt+ (l-i)log(l-i)), (0<f<l, elliptic case), 
a(tlogt-t), (t > 0, flat case), 

a(tlogt- (l + f)log(l + i)), (i>0, hyperbolic case). 

In the elliptic case, u{t) — +oo for t > 1. The Legendre-Fenchel transform of u is 
given by 

'alog(l + e s/a ), (s e E, elliptic case), 
I(s) = { ae s / a , (s G R, flat case), 

-alog(l — e s / a ), (s < 0, hyperbolic case). 

In the hyperbolic case, I(s) = +oo for s > 0. We have i?o = 1 i n the hyperbolic 
case and R Q = +oo in the remaining two cases. The proof is completed by applying 
Theorem [231 □ 
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Proof of Theorem \2.3\ We are going to apply Theorem |2.9| to the polynomial G„(z) 
W n (e l3 n a z). We have / fe> „ = e f3k+aklogn w k . Equation g implies that assump- 
tion (A3) is satisfied with 



u(t) 



a(t\ogt~t), ie[0,l], 
^+00, f > 1. 

The Legendre-Fenchel transform of u is given by 

( ae s / a , s < 0, 



I(s) = 



s > 0. 



Applying Theorem 2.9 we obtain that -/iG„ converges in probability to the required 
limit. A.s. convergence will be demonstrated in Section 5.3 below. □ 



Proof of Theorem{2~b] We apply Theorem [2j| to G n (z) = W(e n a z). We have 
u(t) = a(t\ogt - t) for all t > 0. Hence, I(s) = ae s / a for all s e M. We can 
apply Theorem |2.9| to prove convergence in probability. A.s. convergence will be 
demonstrated in Section [oT4l below. □ 



Proof of Theorem 2.8 Put a — +1 in the case a > 1 and er = —1 in the case a < 1. 
We have u(t) = at a for i > 0. It follows that 



J(r) 

We can apply Theorem |2. 9 




ar > 0, 
or < 0. 



4. Proof of Theorem 12.131 



□ 



4.1. The logarithmic moment condition. The next lemma is the only place in 
our proofs where the logarithmic moment condition appears. 



Lemma 4.1. Let £o>£ij ■ ■ ■ be i.i.d. random variables. Fix e > 0. Then. 
(19) M := sup < +00 a.s. Elog(l + |£ |) < 00. 

fc=0,l,- e 

Proof. For every non-negative random variable X we have 



00 



J^P[X > k] < EX < }_^P[X > k]. 

k=l k=0 

WithA = i log(l+|Co|) it follows that Elog(l+|£o|) < 00 if and only if J2T=i P Mo\ > 
e ek — 1] < 00 for some (equivalently, every) e > 0. The proof is completed by ap- 
plying the Borel-Cantelli lemma. □ 

4.2. Upper bound. Take some z £ D/j o \{0}. Fix an e > 0. We will show that 



(20) 



lim P 

n— yoo 



|G„(z)| >e 



n(/(log| z |)+4e) 



0. 



First, we estimate the tail of the Taylor series (10 1 defining G„. By making e 
smaller we may assume that \z\e 2e < Rq. By assumption (A4) there is A > such 
that for all n > A and all k > An, 

|/ fe , n | < (\z\e^)- k . 
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Lemma |4. 1| implies that for some a.s. finite random variable M, 



(21) 



£,kfk,n'- 



k>An 



< M £ e" k \f Kn \\z\ k < M Y, ^ < M - 



k>An 



k>An 



The last inequality holds if n is sufficiently large. Note in passing that this implies 



that for large n the series (10) converges with probability 1. 



We now consider the beginning of the Taylor series ( 10 ) defining G n . Take some 



5 > 0. By assumption (A3), there is N such that for all n > N and all k < An, 
(22) \fk,n\<(f(k/n) + 5) n . 



It follows from (11) that for alH > 
(23) 



ilog|z|+log/(t)</(log|z|). 



Using (22), (23) and Lemma 4.1 with e/A instead of e we obtain that there is an 
a.s. finite random variable M' such that 



(24) 



0<fe<An 



< 



0<k<An 

< M'e en 



( e $*>*\»\+**m)+6\z\ 



0<k<An 



< M'r 2en 



e— (e I{l °^ +5\z\ J 
< M'e 3e "e" /< - los ' z '- ) 



where the last inequality holds if S — 5(e) is sufficiently small. Combining (21) 
and J24|), we obtain that for large n, 



(25) 



|G„(z)| < M'e"W log l z l) +3£ ) +M. 



Since M and M' are a.s. finite by Lemma 4.1 this implies (20) 



4.3. Lower bound. Fix e > and z £ Dfl o \{0}. We show that 

1 



(26) 



\G n (z)\<e 



n(/(log|*|)-4e) 



= o 



We will use the Kolmogorov-Rogozin inequality in a multidimensional form 
which can be found in [8]. Given a e?-dimensional random vector X define its 
concentration function by 



(27) 



Q{X;r) = sup F[X G D r (a;)], 



r > 0, 



where D r (x) is a d-dimensional ball of radius r centered at x. An easy consequence 



of (27) is that for all independent random vectors X, Y and all r, a > 0, 
(28) Q{X + Y;r) <Q{X;r), Q(aX;r) = Q(X;r/a). 

The next result follows from Corollary 1 on p. 304 of |S]. 
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Theorem 4.2 (Kolmogorov-Rogozin inequality). There is a constant Cd depend- 
ing only on d such that for all independent (not necessarily identically distributed) 
random d-dimensional vectors X±, . . . , X n and for all r > ; we have 



-1/2 



Q{X X 



+ X n ;r)<C d - (j^(l-Q(X fc ;r)) 



\k=l 



The idea of our proof of ( 26 ) is to use the Kolmogorov-Rogozin inequality to show 
that the probability of very strong cancellation among the terms of the series ( 10 ) 
defining G n is small. First, we have to single out those terms of G„ in which the 
coefficient fk, n is large enough. By definition of /, see ( 11 ), there is to € [0, T ] such 
that to log M + log f(to) > /(log \z\) — e. Moreover, by assumptions (A2), we can 
find a closed interval J of length |J| > containing to such that 

/(t)|z|* > e^sI'D-a*, t € J. 

Define a set J n = {k € No : k/n G J}. By assumption (A3) there is N such that 
for all n > N and all k € J n , 

\f k ,n\\z\ k >e n W°^- 3e) . 



-n(/(log|z|)-3e) 



fk . n % 



For k = 0, . . . , n define 
Note that \ak, n \ > 1 for k G J n . Define 

Gn.l = J]] 0-k,n^kj £*n,2 = J]] a k,n^k- 

kej n k<£J n 

By taking the real and imaginary parts we can view the complex random variables 
«fc,n^fc as two-dimensional random vectors. Using (28 1 we arrive at 

(29) P[|G n (z)| <e"( / ( 1 °gl 2 D- 4£ )] <Q(G, M + G„, 2 ;e- £ ") <Q(G„, i; e- £ ™). 

By Theorem |4.2| there is an absolute constant C such that for all r > 0, 

Q(G n tfr)<C'(^(l-Q(a k , n tor))\ < C ■ ( £ (1 - Q(£ fc ; r))\ 

\kej n J \k£j n ) 

Here, the second inequality follows from the fact that |dfc n | > 1. Now, since the 
random variable £o is supposed to be non-degenerate, we can choose r > so small 
that Q(£o; r) < 1. Note that this is the only place in the proof of Theorem 2.9 where 
we use randomness. The rest of the proof is valid for any deterministic sequence 
£oj£i) ■ ■ ■ such that |£ n | = 0(e Sn ) for every S > 0. If n is sufficiently large, then 



< r and hence, 



(30) 



Q(G nA ;e- en ) < Q(G nX ,r) < C^J^ 1 ' 2 < C- 



in 



-1/2 



In the last inequality we have used that the number of elements of J n is larger than 
(| J\/2)n for large n. Taking (29) and (30) together completes the proof of the lower 
bound fl26l). 
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5. Proof of Theorem 12^91 and related results 



5.1. Proof of Theorem 2.9, Recall that /xg„ is the measure counting the zeros 



of G n . Our aim is to show that for every smooth, compactly supported function 

<p : Ofl,, -> R, 

1 



(31) 



tp(z)/iG n (dz) 



p ^ 

n— ¥ oo 



tp(z)fi(dz). 



Let E n be the event G„ = 0. The left-hand side of (31) is not well-defined on 
E n . However, we will argue that lim n _ J . 00 P[-E n ] = 0. By assumptions (Al), (A2) 
and (A3), if n is large, then f k<n > for all k < T a n/2. Also, P[£ = 0] < 1. It 
follows that limn^oo Pf-E^] = 0. In our proof of (31) we may restrict ourselves to 
the complement of E n . 

It is known, see [T2], that on the complement of E n , 

(32) / <p(z)n Gn (dz) = ±- [ A v (z)\og\G n (z)\X(dz). 

Jd Ro ^ Jd Ro 



This is just a restatement of (|15j). From Theorem 2.13 we know that for every z € 
Dr o \{0}, the random variable p n (z) := ^ log |G„(z)| converges to p(z) :— /(log \z\) 
in probability. Assuming for a moment that we can pass to the limit under the sign 
of integral (which will be justified later), we have 

(33) i / V (z)» Gn (dz) A -L / A<p{z)p{z)\{dz). 



This proves Theorem 2.9 but with the formula sp-A(J(log |z|)) for the limiting 



distribution of zeros, as in Remark 2.14 We prove now that the limiting distribution 



of zeros can be given by ([12]) ■ We claim that with /i defined by ( 12 ), 
1 

2r7 



(34) 



Acp(z)p(z)\(dz) 



(p(z)n(dz). 



If / is smooth, then by Green's identity, the left-hand side of ( |34[ ) is equal to 
1 



2rr 



<p(z)A(I(log\z\))\(dz) = 



I"(log\z\) 

v{z) ^u^ x{dz) = 



<p(z)fi(dz). 



Bo 



If / is not smooth, we can find a sequence of non-decreasing, smooth functions 
I' 1: I' 2 t ■ ■ such that I' n < I' and limn-j.oo I' n (s) = I'(s) for all s < logi?o where I' 
is continuous. By dominated convergence, I n (s) :— J_ I' n (t)dt — > Z(s) as n — > oo. 
For each I n we can use the Green's identity as above, and then let n — > oo to 



obtain (34) in full generality. 



5.2. Dominated convergence. It remains to justify the interchanging of the limit 



and the integral when passing from ( 32 ) to ( 33 ) . Recall that a sequence of random 



variables X n is bounded in probability (or tight) if for every e > we can find 
A = A(e) such that P[|A„| > A] < e for all n e N. Recall also that X„ is called 
bounded a.s. if limsup„_ J . oc \X n \ < oo a.s. We need a lemma from [3"I] . 

Lemma 5.1 (Lemma 3.1 in |34j). Let (X,A,v) be a finite measure space. Let 
f n : X — > K be random functions defined on a probability space (Q, £>, P) which are 
jointly measurable with respect to A® B. Assume that for v-a.e. x £ X we have 
f n {x) — >• in probability (reap., a.s.) and that the sequence J x \ f n {x)\ 1+s v{dx) is 



18 



ZAKHAR KABLUCHKO AND DMITRY ZAPOROZHETS 



bounded in probability (resp., a.s.) for some 5 > 0. Then, J x f n {x)v{dx) converges 
in probability (resp., a.s.) to 0. 

Recall that the function if vanishes outside some disk D r , where r < Rq. By 
Lemma 5.1 the passage from (32) to (33 1 is justified if we show that the sequence 
of random variables 

(35) B n := [ \Atp(z)\ 2 \p n {z)-p(z)\ 2 \(dz) 

JjO r 

is tight. Since tp is bounded on D r , we have 

B n <C f \p(z)\ 2 \(dz) + C f \ Pn (z)\ 2 X(dz). 

The first summand on the right-hand side is a finite constant since < /(log \z\) < 
/(log r) . To complete the proof of Theorem 2.9 we need to demonstrate the tightness 
of the sequence 

(36) B„ :=\ ( (log|G„(z)|) 2 A(dz). 



Let us consider first the case in which G„ is a polynomial of degree n and 
/(l) 7^ 0. This applies for example in the setting of Theorem 2.3 Denoting by 
wi n , ■ ■ ■ j w nn the zeros of G n we can write 

G„(z) = Cnfn,n{z ~ W ln ) . . . (z - W nn ). 

Using the inequality of the arithmetic and quadratic means two times we obtain 



2-Trr 2 2 «A, 

B n < —(log |/„.„| + log |£„|) 2 + 



k=l 



log 2 \z - w kn \\(dz). 



Recall that Imin^^ — log = log/(l) is finite. Hence, the first term on the 

right-hand side is bounded in probability. The second term is bounded by a deter- 
ministic constant since J B log 2 \z — w\X(dz) < C for some constant C — C(r) not 
depending on id 6 C. 

In the general setting the proof of tightness of B n is more involved. The main 
difficulty is that log|G„| becomes infinite at zeros of G„. Thus, we have to show 
that with high probability, G„ has not too many zeros. We will prove that B n is 
bounded a.s.: 



(37) 



limsup-B„ < oo a.s. 



If G„(0) = we let r„ be the multiplicity of the zero at 0. Write G* (z) = 
G n (z)/z Tn and define G*(0) ^ by continuity. Let M n (R) = supi z | =R |G* (z)\, 
where R < Rq. First of all, it follows from (25 1 that 



(38) 



Di{R) := limsup — log M n (R) < oo a.s.. 



Note that M' and M there do not depend on argz. Let N n (R) be the number of 
zeros of G* in the disk Dr, counting multiplicities. Denote by a\, n , . . . , ajv„(fl),n the 
zeros of G* in the disk Dr. The Poisson- Jensen formula, see, e.g., [23 Chapter 8], 
states that for every z £ D r and every R > r we have 



(39) 



N n (R) 

log\G* n (z)\ = I n (z;R)+ J2 lo S 

k=l 



R{z - a kiTl ) 



R 2 - a k „z 
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where 
(40) 



I n (z;R) 



2tt Jo 

and Pr is the Poisson kernel: 



log\G* n (Re ie )\P R (\z\,6 - aigz)de 



R 2 - P 2 



R 2 + p 2 - 2Rpcos(p' 



Fix a small 6 > 0. Note that < P r+ s(p,tp) < C for all p < r and ip e [0, 2n], 
where C depends only on r and 6. It follows that I n (z;r + 6) < C log M n (r + 6) for 
all z£D r and hence, by (38), 

(41) 

We will show that 



£>2 := limsup — sup I n (z; r + 6) < oo a.s. 

n— yoo Tl z£D r 



(42) 



limsup —N n {r + 6) < oo a.s. 

n— > oo ^ 



It follows from the Poisson- Jensen formula ( 39 ) with R = r + 26 and z = that 

M„(r + 2<5) 



(43) 



iV n (r + <5) < 



log 



1 

r+28 



lop 



G*(0) 



Indeed, any zero of G* in V> r +s gives contribution at most log to the sum on 
the right-hand side of (39) and 7„(0;i?) < M n (i?). Clearly, G;(0) = f Tn ,nir n - Let 
r = min{n : ^ 0} be the index of the first non-zero £„. Note that r < oo a.s. 
By (Al), (A2), (A3), for sufficiently large n we have f k>n > for all k € [0,T o n/2]. 
It follows that t„ = t for sufficiently large n. Hence, 



(44) 



liminf- log |G;(0)| = a.s. 



Applying ( 38 ) and (44 1 to the right-hand side of ( 43 ) we arrive at (42 1. 
Next we show that 

(45) D3 := liminf — inf I n (z:r + 6) > — oo a.s. 

n->-oo n zSD r 

It is known that I n (0;R) = ^ J Q w log \G^{Re %e )\d6 is non-decreasing in R. Note 
that 7„(0;0) = log|G*(0)|. Hence, by ((44), 



(46) 



liminf — J„(0; R) > —oo a.s. 

n— >oo 71 



From now on we set R = r + 6. Let q n {9) — — log |G* (Re )\ and write 

q+(9) = max(g n (0), 0), = - nrfn(?„(0), 0). 

Then, g„(0) = — %.{^)- Note that there is a constant C > 1 depending only 

on r, (5 such that 1/C < Pr(p, if) < C for all p < r and (/? € [0, 2tt]. We have, for 
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all z S B r , 
2tt 



I n (z;R) 



q+(9)P R (\z\,e-aigz)de- 



2tt 



q-(9)P R (\z\,e -aig z)dS 



> 



> 



1 

CJ 
2ttC 
n 

2irC 
n 



2 71 



2;r 



-C q-{Q)dO 



I n (0;R)- \ C- 



C 



2tt 



I n (0;R)- [C- 



1 \ 2?r 



C7 n 



g+(0)d0 



log M„(B). 



Recalling ( 38 1 and ( 46 ) we arrive at ( 45 ) . 

We are ready to prove that the sequence B n is bounded a.s. Applying the 
inequality of the arithmetic and quadratic means to (39), we obtain 

R(z - a k . n ) 



N n (R) 



(log|G„(z)|) 2 <3(T n log\z\) 2 + 3I*{z;R) + 3N n (R) £ log 2 



k=l 



R 2 - a k ,nZ 



There is a constant C depending only on r, 6 such that for all a € O/j, 



log 2 



R(z - a) 
R? - az 



\{dz) < C. 



Recalling (36 1 we have, for some constant C depending only on r, 5, 

C 



Bn < 



T z n + sup P n (z;R) + N r i(R) 



Recall that r n = r for sufficiently large n. Utilizing (41), (45), (42) we arrive 
at (37). The sequence B n is bounded a.s. and hence, tight. 



5.3. Proof of the a.s. convergence in Theorem |2.3[ Convergence in proba- 
bility has already been established in Section [3j Given n £ N we can find a unique 
jn e N such that j n < n < (j n + l) 3 . Write m„ = j 3 and G n (z) = W„(e^m"z). 
Note that limn^oo m n /n = 1. We will show that -/iG„ converges a.s. to the mea- 
sure with density Q. To this end, we need to prove the a.s. convergence of the 
corresponding potentials. Fix z £ D. We will prove that 



(47) 



1 



log | G„ (z) | 



n— } oo 



Note that G„ satisfies all assumptions of Section 2.5 It follows from ( 25 ) that 



1 



limsup — log |G„(z)| < a\z\ /a a.s. 

n— too tl 



Thus, we have to prove only the lower bound in (47). Fix a small e > 0. It follows 
from (26) and the Borel-Cantelli lemma applied to the subsequence {i 3 }jeN that 
with probability 1 for all but finitely many n £ N, 



(48) 



\G m Jz)\>e r > 
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Let now n be a number not of the form j 3 . We have, by Lemma 4.1 and ^ 



\G n (z) - G m Jz)\ = 



n 

E 

k—m n -\-l 



i k w k e 0k m a n k z k 



<Me 2en e- a(klosk - k) n ak \z\ k . 

k=m n + l 

The function x H ► —a{x log x — x) +ax log n defined for x > attains its maximum, 
which is equal to an, at x = n. Recall that \z\ < 1. Since m n > (1 — e)n and 
M < e en if n is sufficiently large, we have the estimate 

\G n (z) - G m „(z)\ <e^ n e an \z\^ n . 

Since a + log |z| < a^l 1 '", we have, if e > is small enough, 

(49) |G n (z) - G m „(z)| < e (l-^)»(«k| 1/ °-2e) < e m„( a | Z | 1 /°-2s)^ 

Bringing ( [48] ) and ( 49 ) together we obtain that with probability 1 for all but finitely 
many n, \G n {z)\ > e m »H*l /Q - 2£ ). This is the required lower bound in ((47]). 

We are ready to complete the proof. We need to show that ^Mg„ converges a.s. 
to a measure fi with density Q. Take any smooth function ip : C — > R having a 
support contained in D. Write p n {z) = ^ log |G„(z)|. As in (32 1 we have 



n L — ' 



1 

2^ 



A(p(z)p n (z)X(dz). 



We have shown in (47) that for every z € D, p n (z) converges to p(z) = alzl 1 /" a.s. 
Assuming that interchanging the limit and the integral is possible, we arrive at 

1 



(50) 



Sn{f) 



a.s.^ 
n— too 27T 



A(p(z)p(z)\(dz) 



tp(z)u(dz). 



The last step follows from (34). To justify the interchanging of the limit and the 
integral we use Lemma 3.1 of [34j . To use the lemma, we need to show that B n 
defined in (35) (or, equivalently, B n defined by (36)) is bounded a.s. This means 
that limsup,^,^ B n < oo a.s. But we have already verified this in (37). 

Unfortunately, we were able to establish (47) for z £ D only. That is why the 
support of (p was restricted to D. To complete the proof we have only to get rid 
of this assumption. Given a small e > let tp e : C — > [0, 1] be a smooth function 
which is 1 on Di_2e and outside Di_ e . Let N n (r) be the number of zeros of G„ 
inside D r . Then, by 

1 



(51) 



'N n {l-E) > SMs) 



a.s.^ 

n— *-oo 



Let now ip be an arbitrary smooth compactly supported function on C. Write ip = 
fle+^ei where ipi i£ = ip-ip £ is smooth with support in Di_ £ and ipi,,e — ip-(l — ip £ ) 
is smooth with support in C\Oi_2e- Then, by |50|) and (51), 



a.s. 

n—^-oo 



ip llS (z)fi(dz), limsup5„(^ 2 , e )< (l-(l-4e) 1 / Q )||^|| 



Note that ipi >e coincides with ip on Di_2e and hence, J B <pi tS (z)u(dz) converges 
to f n (p(z) u(dz) as £ j 0. Since e > was arbitrary, it follows that S n (ip) = 
S n (ifi, £ ) + S n (ip2, £ ) converges a.s. to L tp(z)[i(dz). This completes the proof. 
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5.4. Proof of the a.s. convergence in Theorem |2.6| Let m n be defined in 
the same way as in the previous proof. Write G n (z) = W(e /3 m"z). We will show 
that the sequence of random measures ^Mg„ converges a.s. to the measure with 
density Q. This implies Theorem 2.6 since lim^oo m„/n = 1. Fix z G C. Wc 
show that 



(52) 



1 



log|G„(s) 



a.s. 

n— ¥ oo 



a\z\ 



l/a 



Note that G„ satisfies the assumptions of Section 2.5 with 7(s) = ae s ' a there. 
From (pBl it follows that 



1 



limsup — log |G„(z)| < a\z\ /a a.s. 

n— too 17> 



We prove the lower bound in (52). Fix small s > 0. It follows from (26) and the 
Borel-Cantelli lemma applied to the subsequence {j }jeti that with probability 1 
for all but finitely many neN, 

|G m „(z)| > e ™"H*l 1/Q - e ). 

However, G rrln (z) — G„(z) by definition. Also, lim n ^oo m n /n = 1. This proves 
the lower bound in (52). The rest of the proof is the same as in Theorem 2.3 but 
we don't need to worry about the case z D. 

5.5. Proof of Theorem 



2.5 



Let W n (z) = X]fc=o £kWkZ k , where Wk is a sequence 
satisfying ^ and ((SI). Assume that Elog(l + |£ |) = oo. Fix e > 0. We will 
show that with probability 1 there exist infinitely many n such that all zer os of 
W n (e~^n~ a z) are located in the disk D £ . We use an idea of Q3]. By Lemma 4.1 
limsupyj^^ ICnl 1 ^" = +oo. Hence, with probability 1 there exist infinitely many n 
such that 



(53) |f n |» > max |£„ 

k — l,...,n — 1 



> max 



3C 



1 1 



Let n be such that (53) holds. By ^ and (53), we have for every z £ C and k < n, 



^n—ksn—k 



i — k 



< C\w n \e pk n ak \^ n 



8 n a 



= C 



For every z such that \z\ > e we obtain 

i—k 



n-l 

E 



(^) 



< c 
1 



n I — 

^ n G^O r ( S 



\k=l 



(3C+1) 



< 



Z 

ePn a 



By (J3j) and (53), the right-hand side of this inequality goes to +00 as n — > 00. For 
sufficiently large n, it is larger that |£o^o|- It follows that for \z\ > e, the term of 
degree n in the polynomial W n (e~ l3 n~ a z) is larger, in the sense of absolute value, 
than the sum of all other terms. Hence, the polynomial ~W n (e~^n~ a z) has no zeros 
outside the disk D F . 
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5.6. Proof of Th eorem |2.10] Start with a measure n satisfying the assumptions 



of Theorem 



2.10 



Define a function / by I(s) — J_ /z(O e r-)dr for s < logi? - 
The integral is hnite by the second assumption of the theorem. Clearly, / is non- 
decreasing, continuous and convex on (— oo, log Ro). For s > logi?o let I(s) = +oo. 
Define /(log Rq) by left continuity. Let now u be defined as the Legendre-Fenchel 
transform of /: 

u{t) = sup(st — I(s)). 

We claim that the random analytic function G n (z) = ^frL n £,kfk,nZ k with ff- t n — 
e -nu(k/n) satis fi es assumptions (Al)-(A4) of Theorem]^ with / = e~ u . By 
the Legendre-Fenchel duality, the function u possesses the following properties. 
Firstly, it is convex and lower-semicontinuous. Secondly, it is finite on the interval 
[0,Tb), where To = limsup t _ i . +00 I(t)/t satisfies To £ (0,+oo]. This holds since / 
is non-decreasing and lim s _ ! ._ co I(s) — by construction. Thirdly, u(t) = +oo for 
t > To and t < 0. This verifies assumption (Al). Fourthly, formula {jTTJ) holds and 
lim t _j. +00 u{t)/t — logi? . This, together with Lemma [4. 1| shows that the conver- 
gence radius of G„ is Rq a.s. and verifies assumption (A4). Finally, u is continuous 
on [0,To) (since it is convex and finite there), and, in the case To < +oo, the func- 
tion u is left continuous at To (follows from the lower-semicontinuity of u). This 
verifies assumption (A2). Assumption (A3) holds trivially with / = e~". 

5.7. Proof of Proposition |2.15| The logarithmic potential generated by the mea- 
sure (J,( K ) is 

U (k) (z):= f log , 1 l f j,W(dw). 
Jc \z-w\ 



We will show that 



'«(«;)- u(0), |z|<e"'(°\ 
(54) U M {z) = - { /(log |^|) +u(k), e M '(°) < \z\ < e u ' (K \ 

K log|z|, \z\>e u ' {K \ 



Here, u'(0) = \im m u' (t) . It follows from d54b that F(z) := U^(z) +/(log|z|) is 



constant and equal to — u(k) on the annulus e" ^ < \z\ < e" ^ which contains 
the support of /i^. For \z\ > e u ^ we have F(z) — I{\og\z\) — K\og\z\ > —u(k) 
by (fill. Finally, for \z\ < e u 'W we have F(z) = /(log \z\) + u(0) - u{n) > -u(k) 



by (11). To summarize, F(z) is constant on the support of fi^ and is at least as 
large as this constant outside the support. By Theorem 3.3 on p. 44 of [24] this 
implies that ^jS k ^ is the equilibrium measure for the logarithmic potential in the 
presence of external field /(log \ z\). The theorem mentioned above is stated in [23] 
for k = 1 only, but it is valid for every k > 0. 



It remains to prove ( 54 ) . The uniform probability distribution A r on the bound- 



ary of D r generates the potential 

(55) / log, — - — 7\ r (dw) = - max{log r, log \z\}. 

JdB r \Z-W\ 

Write F(r) = Z'(logr). Recall that /z^ K ' is a rotationally invariant measure of 
total mass k having support in the annulus e" ^ < \z\ < e u ^ and such that 
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H {K) (B r ) = F(r). Together with ([55) this immediately implies ([54]) for \z\ > e u '^\ 
Suppose now that \z\ < e u Using (55 1 and then integration by parts, we obtain 

u' (k) 



.'(0) 



\og\s\dF{s) = I'(r)r 



«'(0) 



u.'M 



ds 



/'(logs) — =u(k) -it(0). 

u'(0) s 



Here, we used the identity I(u'(t)) = tu'(t) — u(t) for i = and t = k. Suppose 
finally that e u '^ < \z\ < e u '^\ Using Q wc obtain 



U {k) {z) = / \og\z\dF(s) + 



log |fl|dF(«)= /(log |«|) +«(«), 



where the first integral is equal to _F(z)log|z| and integration by parts has been 
used for the second integral. The proof of (54) is completed. 
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